Combined effect of thermal and quantum fluctuations in superconducting 
nanostructures: a path integral approach 



Pedro Ribeiro 

Max Planck Institute for the Physics of Complex Systems, Nothnitzer Str. 38, 01187 Dresden, Germany 

Antonio M. Garci'a-Garci'a 
University of Cambridge, Cavendish Laboratory, JJ Thomson Avenue, Cambridge, CBS OHE, UK 

We study the combined effect of thermal and quantum fluctuations in a zero dimensional super- 
conductor. By using path integral techniques, we obtain novel expressions for the partition function 
and the superconducting order parameter which include both types of fluctuations. Our results 
are valid for any temperature and to leading order in 5/ Aq where S is the mean level spacing and 
Ao is the bulk energy gap. We avoid divergences at low temperatures, previously reported in the 
literature, by identifying and treating non-perturbatively a low-energy collective mode. In the low 
and high temperature limit our results agrees with those from the random phase (RPA) and the 
static path approximation (SPA) respectively. 
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Superconductivity in nano-structures has attracted the 
attention of theorists [l| and experimentahsts Q since 
the early days of the Bardeen-Cooper-Schriffer (BCS) 
theory. Explicit calculations soon showed that, even 
within a mean field formalism, finite size effects had a 
profound impact on the superconducting state. Exper- 
imentally it was also observed 0. 0| that the supercon- 
ducting transition in nanowires and small particles be- 
came broader as the grain size decreases due to ther- 
mal and quantum fluctuations. The use of path inte- 
gral techniques 0| led to a quantitative description of 
thermal fluctuations specially in zero dimensional super- 
conductors where the SPA is applicable. By contrast 
semi-phenomenological models Q that combined quan- 
tum and thermal fluctuations provided only a qualita- 
tive description of the broadening of the transition ob- 
served in superconducting nanowires. The field received 
an important impetus in the mid nineties after the exper- 
iments of Ralph et al. on single, isolated Al nanoparticles 
that showed for the first time that superconductivity 
survived in single particles down to a few nanometers. 
These experiments also stimulated the theoretical inter- 
est in ultrasmall superconductors. At zero temperature, 
the Richardson's formalism [lo| . originally introduced in 
the context of nuclear physics, made possible to find ex- 
act solutions for the low energy excitations of the reduced 
BCS Hamiltonian [ll|. However a theoretical analysis 
that takes into account thermal and quantum fluctua- 
tions simultaneously is still an open problem in the field. 
In [3| this problem was addressed by combining the SPA, 
that models thermal fluctuations, with the RPA, that ac- 
counts for quantum fluctuations to leading order in 6/Aq. 
However it was found that the resulting partition func- 
tion had singularities at low temperature. Progress in 
this problem are specially timely as recent experiments, 
taking advantage of advances in the growth and control 
of nanostructures, have put the basis to test quantita- 



tively the limits of superconductivity in the nanoscale 

H, III, [H, [13, [Hi. 

The main goal of this paper is to put forward a 
theoretical analysis free of divergences and valid at all 
temperatures that combines thermal fluctuations and 
quantum fluctuations to leading order in S/Aq in a zero 
dimension superconductor. These results are of interest 
for other strongly interacting Fermi systems beyond the 
realm of superconductivity in nanostructures. Typical 
examples include hot nuclei (see [Tot and references 
therein), and trapped cold atomic gases Technical 
details of the calculation are postponed to a forthcoming 
publication (23|. Here we summarize the main results, 
their limits of applicability and some technical aspects. 
Moreover we also provide explicit results of the order 
parameter in two simple cases: a constant spectral 
density and a single, highly degenerate, shell. 

Model, approximations and main results 
We consider the BCS Hamiltonian: 

H = '^Eacl^Ca.a (1) 

^ /(a,a')c]j_icLa _ic_a' _iCa',i 



where a, —a label one-particle states related by time re- 
versal symmetry with energies = £-q, <5 is the mean 
level spacing, cr = ±1 is the spin label, g the dimension- 
less coupling constant, I{a,a') = L'^ J d'^r'$-^{r)'$-^,{r) 
with an eigenstate of the one-body problem and L 
the system size. The partition function of the model is 
given by Z = Tr [e~^^] . Ferniionic degree of freedom 
can be integrated exactly by introducing a complex val- 
ued Hubbard-Stratonovich field A(r, r) which results in 
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a partition function, 



Z 



/ 



(2) C{z) = (-.2+4.^)(„,2) 



where is the partition function for free electrons. The 
main goal of the paper is to evaluate ([2|) , including ther- 
mal and quantum fluctuations. The main approxima- 
tions in our calculation are: a) the grain size is zero 
dimensional, namely, the coherence length ^, is smaller 
than the system size. As a consequence, A(t, r) only de- 
pends on imaginary time A(t, r) k, A(t); b) the time 
dependence is sufficiently weak so that an expansion to 
second order is justified. At T = this corresponds 
with the usual RPA around the saddle point solution 
Ag which is valid in the limit 5/Ao < 1; c) the inter- 
acting region in the BCS Hamiltonian is restricted to 
the a narrow interval \~E-d,E-d\ around the chemical 
potential where En is the Debye energy. We assume 
that /(a, a') ^ L'^ J (i'^r^'^(r)\E'^, (r) ^ 1 is independent 
on a, a' and /(a, a') — outside the interactin g re gion. 
This is a good approximation as it was shown in |l4| that 
the a dependence is only important for ^ < L ; d) we as- 
sume that for S/Ao < 1 Coulomb interactions can be 
accounted by a simple redefinition of g. Recent experi- 
ments [l^ suggest that, at least for Pb and Sn grains, this 
is a good approximation up to sizes S ~ Aq ot L ^ 5nm. 

We also note that, for the sake of simplicity, no observ- 
able sensitive to odd-even effects are considered. How- 
ever, even in a grand canonical formalism, it is possible 
to take them into account [2l|, at least for low temper- 
atures, by simply blocking the level closest to the Fermi 
energy. 

The main result of this letter is the following expression 
for valid at all temperatures, that includes simultane- 
ously thermal and quantum fluctuations, 
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? = ysl+e^, g{e) = Y^aH^-^a) is the 
spectral density of the one-body problem, nf{z) = 
,j} -, are the Fermi and Bose function 
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respectively and 
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with r (C) = 2^tanh ) and - J%. For T = 
we recover the RPA results [l6| and, for T ^ Tc, Z is 
given by the SPA of 0. 



Calculation Highlights 

We give an overview of the calculation leading to ([3]) 
with special emphasis on the main differences with 
respect to the techniques of [Ig. A detailed account of 
technical details will be provided elsewhere [23 [. 

The task is to evaluate simultaneously the contribution 
to the partition function Z of thermal fluctuations, taken 
into account by integrating exactly over the static compo- 
nent of A(t) (SPA) and quantum fluctuations, arising as 
small (imaginary) time dependent Gaussian corrections 
(RPA) to SPA. Previous approaches to this problem (l6| 
considered indeed small corrections to a static solution 
A(0), A (r) A(0) -I- SA (r) where SA (r) is the second 
derivative matrix for the fluctuations 6 A (r). Then A(0) 
is integrated out exactly but the integral over SA (r) is 
carried out in the Gaussian approximation only. It is 
therefore assumed that any small correction around any 
A(0) is still a local minimum of the action, namely, the 
real part of the eigenvalues of 6 A (t) is always positive. 
However it was found in [16| that some eigenvalues of 
6 A (t) acquire a negative real part as the temperature 
is lowered. As a consequence divergences occur and the 
theory breaks down, preventing thus the combined study 
of quantum and thermal fluctuations. Divergences in this 
context usually suggest the existence of a collective zero 
mode that must be treated non-perturbately. 

In order to identify this collective mode we separate 
phase and amplitude fluctuations by using polar coor- 
dinates A(t) = s(T)e*'^(^), with s(t) = so + Ss{t) and 
(4) 4'i''') ~ ooT + ^0 + '5(/)(t), where (5(/)(t) and (5s(t) are 
small fluctuations around the static values so , (/>o and 
ao = (M € Z) accounts for phase conflgurations 

with non-trivial winding numbers. Only the ao = is 
considered, contributions from different values of ao, 
known to be related to odd-even effects [2l| are not 
addressed here. The identiflcation of this collective 
mode is nevertheless crucial, if treated perturbatively 
it will lead to the negative eigenvalues and divergences 
observed in [l^. By contrast, following the above 
decomposition, the eigenvalues of (50(t) and 8s(t) in 
our case have always a positive real part and therefore 
no divergences arise. We can then treat separately the 
collective mode and integrate exactly over the static 
phase 00 . This is the key difference between our method 



and that of [l^. As a consequence our results provide 
a quantitative, free of divergences, description of the 
combined quantum and thermal fluctuations at any 
temperature. 

Results 

The natural order parameter for the superconducting 
transition is the connected pair correlation function 

= (5'^)^Eaa' (cL'ici'-iCa-iCai)^. An explicit 
expression for Ac is obtained in a standard way by 
adding source terms to the action ([2]) and deriving with 
respect to them, 



(6) 
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Figure 1: Ac(r), ([6]), for a constant spectral density = 
1/5. Ac(r) combines thermal and quantum fluctuations. It 
reduces to the RPA (SPA) for T -^TciT ■> Tc) . 
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defined as 
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and the average ((...)) is 
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In the literature other parameters have been 
considered to study deviations _from mean- field 
results: for example ((*o/ 



(gsy 



The latter can be simply 
A2 - gS {5g) deg (e) 
For simplicity we assume Ap - 



and 



A2 



A^ 



related 
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to ^ by 

{0'))-nf [ef]. 

A as other terms in © 
do not play a significant role and make the calculation 
slightly more involved. Ac becomes the bulk gap for 
(5 — ?> 0, and it is expected to be closely related to the 
spectral gap at finite 5. We focus on two specially simple 
situations: a) a constant spectral density, b) only one 
level, usually called shell, in the interacting region 
with a degeneracy Ni ^ 1 such that i5/Ao <C 1 where 
5 = 2E]j /Ni . Physically this corresponds to a spherical 
or cubic grain in which, due to geometrical symmetries, 
the spectrum is highly degenerate. Other geometries can 
be easily studied but calculations are more involved. We 
postpone this study to a future publication (23| . 



Constant spectral density 
In this case g{e) = 1/6 and the partition function ^ 
cannot be simplified further so we carry out the calcu- 
lation of Ac (El numerically. In FigUwe depict Ac(T) 
for different values of 6. As was expected, no divergences 



arise at low temperatures. For zero temperature Ac(0) 
is equal to the RPA result [l^ that predicts a leading 
correction Ac(0) = Ao(l -I- aS/En) with a a constant of 
order the unity. For T :^ Tc Ac agrees with the SPA 01 
that describes thermal but not quantum fluctuatio ns ( see 
Fig. [2]). Results from Richardson's formalism [lol . 22 1 at 
T ~ are similar but a direct comparison is not possible 
as Ac is not exactly the spectral gap. In Fig. [2]we depict 
the difference between ([7]) and the SPA prediction. De- 
viations at low temperatures are mostly due to the RPA 
correction, however it is clearly observed that, for inter- 
mediate temperatures, differences from SPA results in- 
crease as a consequence of the combined effect of thermal 
and quantum fluctuations. Previously this region was not 
accessible to analytical calculations. We note that the 
observed enhancement of Ac by quantum and thermal 
fluctuations is not an indication that superconductivity 
is more robust. In fact fluctuations always weaken long 
range order causing phase slips and the broadening of 
the transition. The gap is enhanced because fluctuations 
induce pairing in circumstances which are not allowed by 
a mean field formalism. 

Shell models 

The calculation of the partition function greatly simpli- 
fies by assuming that there are only two degenerate levels 
(shells) in the interaction region. We note that quantum 
fluctuations are still small, and therefore our formalism is 
still applicable, provided that the degeneracy of the level 
Ni/2 is large enough such that 5 ^ Aq where in this 
case S = 2Eo/Ni. With this simpliflcation it is possible 
to flnd an explicit expression for Ai . For two shells with 
energy at ieo (*•£. g{e) = ^[5 {e - Eq) + 5 {e + Eq)]) 
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Figure 2: Difference between Ac{T) ^ and the SPA pre- 
diction Q which only takes into account thermal fluctuations 
for different S's. It is assumed that p(e) « 1/S. For low 
temperatures the difference is just the usual RPA correction 
that describes quantum fluctuations at T = 0. However the 
peak observed close to Tc is due to the non-trivial interplay of 
thermal and quantum fluctuations which is beyond the reach 
of SPA and RPA separately. This is the flrst time that this 
region is accessible to analytical techniques. 



^o,-4i in ^ are given by, 
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For r = the first correction to the mean-field 
result coincides with the RPA prediction, Ac = 

2" 




1+ 



where 



In the limit T > Tc , it is also 
possible to obtain explicit expressions of Ac by expand- 
ing the action in powers of sq. To the lowest order in S, 



the SPA result A 



c 



5gtanh2(^)coth2(^) 
/3[l-tanh(^) coth( 



ered. Higher order terms include deviations from SPA 
due to quantum fluctuations. 

Natural extensions of this work include the calculation 
of thermodynamical observables, odd-even effects, mag- 
netic susceptibilities and the differential conductance, the 
outcome of STM experiments. It would also be interest- 
ing to study in detail a hemispherical particle as this ge- 
ometry can be investigated experimentally (l3| . For the 
sake of clearness we have focused here in developing the 
formalism. Some of these applications will be discussed 
in future publications [23 1. 

We have studied the combined effect of thermal and 



quantum fluctuations in a zero dimensional supercon- 
ductor. For the first time we have obtained explicit ex- 
pressions for Z and Ac(r) valid for all temperatures 
and to leading order in S/Aq. For intermediate tem- 
peratures both fluctuations contribute substantially to 
Ac{T). These results provide a solid theoretical frame- 
work to describe quantitatively pairing in confined ge- 
ometries at finite temperature beyond the mean field ap- 
proximation. A problem of current interest in condensed 
matter, nuclear and cold atom physics. 
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